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INTRODUCTION 


In  this  report  we  describe  a  model  to  predict  the  pressure  field  developed  in  shallow  water  by  a  linear  array 
of  equally  spaced,  discrete  explosive  charges  placed  near  a  sand  bottom.  The  method,  suggested  to  the  author  by 
Gaspin,*  has  a  theoretical  foimdation  that  rests  on  the  works  of  Arons  and  Yennie,^  Rosenbaum,^  Cagniard  et  al.,^ 
Snay,^’^  and  Britt,^’^’*  who  treated  in  considerable  detail  the  pressure  field  produced  by  a  single  underwater  charge 
in  the  presence  of  free  surface  and  bottom  interfaces. 

Presently,  the  pressure  field  model  is  an  empirical  model  based  upon  data  from  the  Eglin  Air  Force  Base  test 
site.  This  enables  the  model  to  be  used  to  reliably  predict  the  pressures  associated  with  mine  vulnerability  tests 
conducted  at  Eglin,*  but  die  model  has  not  been  validated  for  other  venues.  The  model  is  capable  of  extension  to  other 
locations  by  making  use  of  the  physical  theory  cited  above,  principally  that  developed  by  Britt  and  Snay.  The 
necessary  physics,  however,  has  not  been  built  into  the  current,  interim  model. 

The  intention  of  this  modeling  effort  was  to  produce  a  computationally  efficient  code  diat  would  describe  the 
pressure  histories  at  arbitrary  field  locations  to  levels  of  accuracy  adequate  for  nondeterministically  modeling  Eglin 
mine  response  data.  The  quasi-impulsive  nature  of  the  mine  kill  response  mechanism,  which  integrates  the  finer  details 
of  the  pressure  pulse  signature,  and  the  presence  of  random  scatter  in  the  mine  kill  data  suggested  that  an  adequate 
mine  kill  prediction  model  could  be  fashioned  using  approximate  classical  modeling  techniques.  The  primary 
modeling  objective  was  to  predict  mine  kill  probability  as  a  function  of  variables  that  were  systematically  controlled 
in  the  tests  with  enough  of  the  underlying  physics  included  to  permit  extrapolation  and  interpolation  of  Ae  model  to 
untested  conditions. 

Given  the  caveat  of  site  dependency,  the  model  presented  below  will  be  shown  to  possess  "accuracy"  by  various 
measures  that  is  perhaps  surprising  in  view  of  the  model  simplicity.  Comparisons  are  made  with  gage  pressure 
histories  at  all  interesting  in-water  gage  positions  used  in  the  1995  Eglin  tests.  Eglin  pressure  histories  obtained  at 
gage  positions  2  inches  above  the  sand-water  interface  are  used  to  calibrate  the  model  and  to  confirm  the  basic 
physical  process  by  which  line  charge  pressures  arise. 
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MODEL  DEVELOPMENT 


Let  the  origin  of  a  right-hand  coordinate  system  be  on  the  water  surface  with  the  z  axis  positive  downward  and 
the  jc  and  y  axes  as  illustrated  in  Figure  1 .  We  consider  a  linear  array  of  bulk  charges  lying  in  the  y  =  0  plane  parallel 
to  the  X  axis  at  a  depth  dg.  Let  the  first  charge  be  located  at  x  =  0  with  the  remaining  N  charges  extending  in  the 
positive  X  direction  at  spacings  of  Ax^j^g.  Thus,  for  /  ;=  1 . .  .N,  the  ith  charge  has  coordinates 

Xc,  :=  0  -  l)-^chg»  Vc,  :=  (1) 


Assuming  detonation  of  the  charge  at  the  origin  first,  the  detonation  time  of  the  ith  charge  is 


^det,  •“  ^0 


(i  -  l)  Ax, 


det 


'"dct 


(2) 


where  is  the  temporal  origin,  here  taken  as  zero,  Ax^jgj  is  the  length  of  the  detonating  cord  between  charges,  and 
is  the  speed  of  the  detonation  wave  along  the  detonating  cord. 


Figure  1.  Line  Charge  Geometry 
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Let  a  pressure  gage  be  located  at  a  point  with  coordinates  X,  Y,  Z.  As  described,  e.g.,  by  Britt  and  Snay,® 
we  will  assume  that  the  gage  signal  originating  from  each  charge  of  the  line  charge  array  consists  of  three  additive 
components:  (1)  a  direct  component,  (2)  a  surface-reflected  component,  and  (3)  a  bottom-reflected  component. 
The  two  reflected  components  will  be  modeled  using  image  charges  as  shown  in  Figure  2.  We  will  assume  that 
the  complete  pressure  history  produced  by  the  line  charge  array  at  this  gage  position  can  be  represented  by 
additively  combining,  at  each  time  of  interest,  these  three  components  for  all  charges  of  the  array. 


Figure  2.  Ray  Geometries 


The  location  of  the  surface-reflected  image  of  the  ith  charge  is  given  by  coordinates 

X  :=  x^,  y,  :=  y.,  and  z,  :=  z  , 

Si  •'Si  •'  Cf  Sf  Ci 

while  those  of  the  bottom-reflected  image  are  given  by 

:=  Yb,  ■=  Yc;  and  h,  ■= 


Consequently,  the  direct  line-of-sight  slant  range  from  the  ith  charge  to  the  gage  is 


'•d,  +  i¥-y^/  -  (Z-z^/ 


(3) 


(4) 


(5) 
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that  from  the  surface-reflected  image  of  the  ith  charge  is 


+  iY-y,/  +  iZ-z,/  , 


and  that  from  the  bottom-reflected  image  of  the  ith  charge  is 


(6) 


V 

(7) 


Thus,  the  times  of  arrival  at  the  gage  of  the  direct,  surface-reflected,  and  bottom-reflected  signals  from  the  ith 
charge  are  given  by 


and 


(8) 


(9) 


h  ^det  > 


(10) 


respectively,  where  c  is  the  velocity  of  the  shock  waves.  We  shall  use  the  acoustic  approximation  and  assume  c 
is  equal  to  the  water  sound  speed.  For  this  we  employ  a  formula  published  Snay  and  Kriebel^  and  calculate  c  using 
the  program 


r-|-(Temp,-32) 

^-s»v 

1449.3  +r(4.57191  -  r(0.044524  +  7(0.0002604 -  70.0000079851)))... 
+(5-  35)(1.39807  +((5-  35)0.00169202-  7(0.011244  +  70.00000077872))) 

- 

^  0.3048 


(11) 


¥ 


which  converts  their  formula  from  metric  to  English  units  and  provides  c  in  feet  per  second  given  water 
temperature  in  degrees  Fahrenheit  and  salinity  in  parts  per  thousand. 

We  shall  assume  the  direct,  bottom-reflected,  and  surface-reflected  shock  waves  from  all  charges  to  be 
spherical  with  exponential  time  dependence  and  overpressures  of  the  form 

‘-t, 

P(r,t)  «  i»^(r)e‘  ® 


r 
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where  P^r)  is  the  shock  front  overpressure  at  a  distance  r  from  the  charge,  6  is  the  decay  constant,  and  is  tiie 
time  of  arrival  at  the  gage.  For  directly  arriving  shocks  the  proportionality  coefficient  is  taken  to  be  unity.  For 
reflected  shocks,  die  proportionality  coefficients  will  depend  upon  characteristics  of  the  transporting  and  reflecting 
media. 

The  theory  of  surface-reflected  waves  (see  Snay  and  Kriebel^  indicates  that  compression  waves  will  become 
tension  waves  upon  reflection  from  a  free  surface.  Hence,  the  proportionality  coefficient  for  surface-reflected 
waves  will  be  negative  in  sign  so  that  the  amplitude  of  the  reflected  wave  will  be  subtracted  from  the  prevailing 
pressure  that  exists  prior  to  the  arrival  of  the  surface-reflected  pulse.  Since  water  supports  little  pure  tension, 
however,  we  shall  not  allow  the  absolute  pressure  to  fall  below  zero.  In  this  interim  model  we  employ  a 
proportionality  coefficient  of  -1  for  all  surface-reflected  waves.  The  theory  indicates  that  this  is  an 
oversimplification,  particularly  for  waves  reflected  at  shallow  angles.  Eglin  pressure  measurements  show,  however, 
that  gage  pressures  are  strongly  reduced  upon  arrival  of  the  surface-reflected  waves  so  that  their  strengths  may 
be  unimportant.  Experience  with  the  code  suggests  good  agreement  with  pressure  records  is  obtained  with  the  -1 
coefficient  for  gage  positions  near  the  sand-water  interface.  A  more  refined  treatment  of  the  surface-reflected 
waves  may  be  necessary  at  other  gage  positions  and  for  shallower  water. 

The  theory  of  bottom-reflection  is  arguably  more  complicated  than  that  for  free  surfaces.  Britt  has  achieved 
remarkable  success,  however,  in  modeling  shock  wave  reflections  from  complex,  multi-layered  bottom  materials. 
His  REFMS  code  relies  principally  on  linear  shock  propagation  theory,  but  recent  enhancements  of  the  code  have 
extended  it  to  highly  nonlinear  applications.*  Britt  and  Snay  have  suggested  that  reflections  from  materials  that 
lack  shear  strength,  such  as  sand,  can  be  modeled  using  the  simpler  theory  for  non-rigid  or  "liquid"  bottoms.^  The 
reflections  of  plane  front  shock  waves  from  plane  surfaces  are  often  described  by  relationships  like  that  of 
Equation  12  with  coefficients  of  proportionality  that  are  referred  to  as  "reflection  coefficients."  We  shall  employ 
such  a  term  in  the  present  model  and  denote  it  as  4>.  The  reflection  coefficient  $  can  range  in  value  from  that 
appropriate  for  a  fiee  surface  (-1)  to  that  appropriate  to  a  rigid  surface  (+1).  The  theory  shows  that  $  is  a  function 
of  the  reflection  angle,  the  densities  and  wave  speeds  of  the  media,  and  also  of  time.  Negative  values  may  occur 
at  small  angles  of  reflection,  even  \niien  the  density  and  sound  speed  of  the  bottom  material  are  greater  than  those 
of  water.  In  this  interim  model  we  choose  $  empirically  and  represent  it  as  a  function  of  the  reflection  angle. 

Using  simple  free  field  shock  wave  "similitude"  equations  (see,  e.g.,  Cole,^  p.  235),  we  calculate  shock  front 
pressures  for  the  ith  charge  at  ranges  corresponding  to  the  direct,  surface-reflected,  and  bottom-reflected  waves 
as 


respectively,  where  Kp  and  ap  are  empirically  determined  coefficients  that  depend  upon  explosive  type  and  W  is 
the  weight  of  the  ith  explosive  charge.  Similarly,  time  constants  for  the  three  component  waves  are  represented 
«  by 
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Under  the  assumptions  of  the  model,  the  pressure  history  at  the  gage  at  time  t  can  be  represented  by  the 
following  function  which,  except  for  A  and  n,  uses  variables  defined  above.  The  quantity  i4  is  an  ordered  n  x  3 
matrix,  each  row  of  which  consists  of  a  wave  arrival  time  (column  1),  the  index  of  the  charge  from  which  the 
wave  originated  (column  0),  and  an  index  value  equal  to  0,  1,  or  2  that  identifies  the  row  arrival  time  as  a  direct, 
bottom-reflected,  or  surface-reflected  arrival  (column  2)  respectively.  The  rows  of  A  are  placed  in  order  of 
ascending  arrival  times.  Usually  n  is  set  to  a  value  less  than  3N.  To  reduce  calculation  times,  a  truncation  function 
is  used  to  eliminate  those  rows  in  A  after  a  specified  number  of  surface-reflection  arrivals.  For  example,  in  the 
figures  shown  in  the  next  section  the  calculations  were  stopped  after  the  arrival  of  20  surface  reflection  pulses. 


P(t,  n.  A,  Pyg,  Pyj,  0^,  0j,  0^,  tj,  := 


sum- 14 
for  j€\...n-  1 


sum  -  sum  +  P^-e  *"*  if  [(^  =  o]- s: 


|sum-sum  +  ^-Pyj^e  if 

sum-sum  -  Pyj^-e  if  [(,4^2^)y  = 

siun-O  if  sum<0 
sum 


(15) 
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DISCUSSION  OF  MODEL  PERFORMANCE 


Figure  3  shows  a  pressure-time  curve  calculated  using  the  model  for  an  M58  line  charge  array  as  it  was  configured 
in  the  1995  Eglin  mine  vulnerabihty  tests.  This  linear  array  was  composed  of  198  2.5-pound  C-4  charges  spaced  at  6-inch 
intervals.  The  center  of  each  charge  of  the  array  was  located  1.5  inches  off  the  bottom  in  water  10  feet  deep.  The  Jf;  Y,  and 
Z  coordinates  used  in  the  calculation  placed  the  gage  along  the  perpendicular  bisector  of  the  array,  2  inches  from  the 
bottom,  and  at  a  standoff  of  18  feet.  For  simplicity,  a  constant  reflection  coefficient  value,  $  =  -0.996,  was  used  rather 
than  some  more  complicated  ftinction  of  reflection  angle.  This  “constant  <I>  model”  is  compared  with  a  “variable  <1>  model” 
later  in  this  report. 

The  most  prominent  features  of  the  calculated  pressure-time  curve  are  the  pressure  spikes  and  ramp  buildup.  Both 
features  are  due  to  the  interplay  between  direct  and  bottom-reflected  waves.  The  plus  symbols  shown  along  the  top  of  the 
display  are,  from  bottom  to  top,  the  arrival  times  of  bottom-reflected  shocks  (brown),  direct  wave  arrival  times  (blue), 
and  the  surface-reflected  arrival  times  (red)  respectively.  Each  pressure  spike  is.  caused  by  the  arrival  of  a  direct  pulse 
followed  almost  immediately  by  a  bottom-reflected  tensile  wave.  As  we  shall  see,  the  spikes  have  very  short  durations  of 
about  half  a  mia-osecond.  Because  of  the  value,  about  99.6%  of  the  direct  wave  pressure  amplitude  is  removed  by  the 
arrival  of  the  bottom  reflection.  The  ramp  buildup  is  caused  by  the  accumulation  of  the  0.4%  residuals  as  more  subtracting 
pairs  reach  the  gage.  Finally,  the  ramp  is  eliminated  by  the  arrival  of  surface-reflected  waves.  Because  of  the  short  duration 
of  the  spikes,  the  ramp  pressure  is  the  primary  source  of  the  loading  pressure  on  a  target. 


Figure  3.  Calculated  Pressure-Time  Curve 
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Figured.  Leading  Pulses  Figured.  Comparison  of  Averaged  Pulses 

with  Measurements 


To  obtain  the  red  curve  of  Figure  6,  the  pulses  on  the  surface  of  the  gage  at  each  instant  of  time  were  numerically 
averaged.  We  assumed  an  “effective”  gage  diameter  50%  larger  than  the  actual  diameter.  The  need  for  such  a  correction 
is  well-known  in  the  measurements  community.  This  resulted  in  an  approximately  90%  reduction  in  the  amplitude  of  the 
short  duration  spikes.  While  the  amplitude  of  the  leading  pulse  does  not  match  the  data,  the  spread  of  the  leading  pulse 
is  in  good  agreement  with  measurements.  Because  the  amplitudes  of  the  calculated  and  measured  pressures  at  later  times 
are  in  excellent  agreement  it  appears  reasonable  to  speculate  that  spikes  such  as  those  calculated  may  actually  exist,  but 
are  not  measurable  in  detail  with  the  gages  employed. 

In  their  report  on  free  surface  reflections  Snay  and  Kriebel  suggested  that  for  shock  wave  pressures  of  less  than 
10,000  psi  and  ray  reflection  angles  of  not  less  than  30  degrees,  use  of  the  acoustic  approximation  should  result  in  a  vahd 
representation  of  the  physics.^  The  reflection  angles  encountered  in  shallow  water  for  targets  and  line  charges  that  rest 
on  the  bottom  are  much  less  than  30  degrees;  hence,  one  might  expect  the  approximation  to  fail  under  these  conditions. 
Yet  quite  reasonable  results  have  been  obtained  with  the  present  model.  Figures  7  and  8  illustrate  wave  front  geometries 
6  ms  after  detonation  of  the  line  charge  systems  of  current  interest.  These  depict  graphically  how  the  waves  of  particular 
charges  run  within  the  pressure  fields  of  adjacent  charges.  For  the  M58  and  continuous  Sabre  systems,  shock  front 
pressures  drop  below  10,000  psi  at  standoffs  of  about  3  feet,  and  for  the  larger  discrete  Sabre  charges,  the  crossover 
occurs  at  about  4.7  feet.  These  distances,  shown  by  dashed  red  lines  in  Figures  7  and  8,  indicate  that  the  acoustic  wave 
speed  assumption  is  a  reasonable  approximation  ^ou^out  most  of  the  region  of  interest  defined  more  clearly  below. 
This  is  also  borne  out  by  the  agreement  between  the  times  of  surface  cutoff  calculated  and  observed.  Where  the  difference 
between  shock  and  sound  speeds  may  have  significant  effect  is  in  regions  where  arrivals  are  densely  packed,  such  as  near 
the  leading  edge  of  the  pulse  train.  It  is  suspected  that  the  failure  qf  the  model  to  predict  the  initial  peak,  as  well  as  other 
wave  peak  phenomena,  may  be  due  to  such  an  effect. 

In  Figure  9  we  compare  the  specific  impulse  of  the  calculated  pressure  pulse  with  that  determined  from  the  measured 
pressure-time  data  (again  shifted  by  65  ps).  Integration  of  the  calculated  pulse  train  includes  the  spikes,  although  because 
of  the  narrowness  of  the  spike  widths,  the  dominant  contribution  to  the  impulse  derives  from  the  ramp  pressures.  The  close 
agreement  suggests  that  the  model  would  do  well  in  predicting  the  response  of  quasi-impulsively  driven  systems.  Details 
of  the  pressure  history,  such  as  the  locations  of  spikes  and  other  temporally  local  fluctuations,  are  of  course  less  apparent 
in  the  integral  curve. 
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Figure  7.  D-Sabre  Wave  Fronts 


4d 


Figure  8.  MSS  and  C-Sabre  Wave  Fronts 


TIME  (S) 


Figure  9.  Calcuiated  Impulse-Time  Curve 
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MODEL  CALIBRATION 


In  the  above  calculations,  a  constant  reflection  coefficient  value  was  employed  to  calculate  tiie  strengtiis  of  the 
bottom-reflected  waves  from  all  charges  of  the  anay.  Theory  predicts,  however,  that  $  is  not  a  simple  constant,  but 
a  dependent  variable  that  changes  with  time,  the  angle  of  reflection,  and  other  conditions  that  may  change  from  charge 
to  charge,  gage  to  gage,  and  one  problem  scenario  to  another.  In  a  truly  predictive  model  a  functional  representation 
of  $  that  is  based  on  the  underlying  physics  and  experimental  data  is  required.  It  is  needed,  in  particular,  for 
applications  to  mine  vulnerability  because  the  loading  pressure  generated  by  a  line  charge  is  a  sensitive  function  of 
bottom  standoff.  Because  the  theory  is  mathematically  difficult  and  not  wholly  developed  for  these  problems,  we 
undertook  an  empirical  study  as  an  initial  step  in  defining  a  reflection  coefficient  function  for  use  in  a  variable  O 
model. 

The  success  shown  above  of  the  constant  model  in  representing  the  pressure-time  and  impulse-time  curves 
in  the  M58  case,  which  was  the  first  case  studied,  prompted  us  to  apply  die  constant  O  model  to  other  gage  positions 
and  to  other  shots.  It  was  of  interest  to  examine  how  die  $  values  that  were  required  to  match  the  impulse 
measurements  varied  from  one  set  of  conditions  to  another.  Calculations  were  performed  for  all  pressure  gage 
positions  2  inches  from  the  bottom  in  the  Eglin  1995  tests  as  well  as  for  some  of  the  gage  placement  that  were  closer 
to  die  free  surface.  In  the  1995  Eglin  tests  each  of  the  three  line  charge  systems  of  interest  (M58,  D-Sabre,  and 
C-Sabre)  had  been  fired  in  water  both  3  and  10  feet  deep  to  provide  a  total  of  sbc  shots.  In  each  shot  the  total  number 
of  gage  positions  placed  2  inches  above  the  sand  was  eight,  and  four  gages  were  placed  at  shallower  positions  that 
depended  upon  the  water  depth.  Table  I  lists  the  conditions  of  the  1995  Eglin  tests  feat  are  analyzed  in  this  report  and 
also  serves  as  a  guide  to  figures  showing  comparisons  between  pressure  calculations  and  test  data. 


Table  /.  Conditions  and  Figure  Numbers  of  Analyzed  1995  Eglin  Tests 


System 

Shot  No. 

Water  depth 

(ft) 

Standoff 

(ft) 

Gage  height 

Figure  Nos. 

M58 

4283 

3 

6 

2  in 

13a,  A-la 

M58 

4283 

3 

8 

2  in 

13b,  A-lb 

M58 

4283 

3 

10 

2  in 

13c,  A-lc 

M58 

4283 

3 

13 

2  in 

13d.  A-ld 

M58 

4286 

10 

12 

2  in 

10a,  A-4a 

M58 

4286 

10 

15 

2  in 

10b,  A-4b 

M58 

4286 

10 

18 

2  in 

10c,  A-4c 

M58 

4286 

10 

18 

6ft 

lOg 

M58 

4286 

10 

22 

2  in 

lOd,  A-4d 

M58 

4286 

10 

26 

6fl 

lOh 

M58 

4286 

10 

30 

2  in 

lOe,  A-4e 

M58 

4286 

10 

35 

2  in 

lOf,  A-4f 

D-Sabre 

4284 

3 

6 

2  in 

14a,  A-2a 

D’Sabre 

4284 

3 

8 

2  in 

14b,  A-2b 

D-Sabre 

4284 

3 

10 

2  in 

14c,  A-2c 

D-Sabre 

4284 

3 

13 

2  in 

14d,  A-2d 

D-Sabre 

4284 

3 

18 

2  in 

14e 

D-Sabre 

4284 

3 

23 

2  in 

14f.  A-2e 

D-Sabre 

4284 

10 

12 

2  in 

11a,  A-5a 

D-Sabre 

4284 

10 

15 

2  in 

11b,  A-5b 

D-Sabre 

4284 

10 

18 

2  in 

lie,  A-5c 

11 


IHTR  2058 


Table  I— Continued 


System 

Shot  No. 

Water  depth 

(ft) 

standoff 

(ft) 

Gage  height 

Figure  Nos. 

D-Sabre 

4284 

10 

22 

lid,  A-5d 

D-Sabre 

4284 

10 

30 

He.  A-5e 

D-Sabre 

4284 

10 

35 

2  in 

Ilf,  A-5f 

C-Sabre 

4285 

3 

6 

2  In 

15a,  A-3a 

C-Sabre 

4285 

3 

8 

2  in 

15b,  A-3b 

C-Sabre 

4285 

3 

10 

2  in 

15c,  A-3c 

C-Sabre 

4285 

3 

13 

2  in 

A-3d 

C’Sabre 

4285 

3 

23 

2  in 

15e.  A-3e 

C-Sabre 

4285 

3 

28 

2  in 

15f.  A-3f 

C-Sabre 

4288 

10 

12 

2  in 

12a,  A-6a 

C-Sabre 

4288 

10 

15 

2  in 

12b.  A-6b 

C-Sabre 

4288 

10 

18 

2  in 

12c,  A-6c 

C-Sabre 

4288 

10 

22 

2  in 

12d.  A-6d 

C-Sabre 

4288 

10 

30 

2  in 

12e,  A-6e 

C-Sabre 

4288 

10 

35 

2  in 

12f,  A-6f 

Comparison  of  Constant  <t>  Model  Calculations  with  Measurements  in  Water  10  Feet  Deep 

The  10-foot  water  depth  results  are  shown  in  Figures  10  through  12.  In  Figure  10  calculated  quantities 
(shown  in  red)  are  compared  with  measured  quantities  (shown  in  blue)  for  shot  4286  of  the  M58  system.  The 
results  that  jqjpeared  above  in  Figures  3  and  4  are  shown  again  in  Figure  10c.  Similar  comparisons  for  shot  4287 
of  the  discrete  Sabre  line  charge  are  shown  in  Figure  1 1  and  for  shot  4288  of  the  continuous  Sabre  system  in 
Figure  12. 

Each  figure  title  specifies  a  gage  letter  that  identifies  the  gage  cable  and  gage  position.  In  all  10-foot  water 
depth  cases  the  cable  A  gage  was  at  a  horizontal  standoff  of  12  feet  from  the  line  charge  axis,  the  B  gage  was 
at  15  feet,  the  C  and  D  gages  were  at  18  feet,  the  E  gage  was  at  22  feet,  the  F  and  G  gages  were  at  26  feet,  the 
H  and  I  gages  were  at  30  feet,  and  the  J  gage  was  at  35  feet.  All  gages  were  located  2  inches  from  the  sand 
bottom  in  the  plane  perpendicular  to  and  bisecting  the  line  charge,  with  the  exception  of  the  D,  G,  and  I  gages, 
which  were  in  the  perpendicular  plane  and  6  feet  from  the  sand  bottom  (4  feet  from  the  free  surface).  Note  that 
the  ordering  of  the  figures  below,  within  each  shot  group,  is  not  always  alphabetical.  The  records  of  shallow-depth 
gages  and  of  faulty  gages  are  usually  omitted. 

Figures  10a  through  lOf  indicate  that  a  constant  value  of  can  be  found  that  will  produce  good  to  excellent 
agreement  between  the  model  impulse  calculations  and  the  impulse  curves  derived  from  the  experimental  pressure 
measurements  at  the  near-bottom  gage  positions  of  shot  4286.  The  model  also  reliably  reproduces  the  ramp 
pressure  history  and  ramp  cutoff  time,  which  coincides  with  the  first  arrival  of  a  free  surface-reflected  pulse. 
Inability  to  predict  the  leading  spike  is  particularly  apparent  in  Figure  lOf,  where  pressures  are  low  and  the  leading 
spike  causes  a  significant  jump  in  the  initial  impulse.  At  such  low  pressures,  however,  gage  drift  may  also  be  a 
contributing  factor  in  the  disagreement.  As  suggested  earlier,  it  is  likely  that  the  initial  spike  is,  in  part,  a 
consequence  of  the  dense  packing  of  arrivals  at  the  leading  edge  of  the  pulse  train.  These,  when  integrated  over 
the  area  of  the  gage  crystal,  may  create  the  illusion  of  a  broad  leading  spike.  It  may  also  involve  the  buildup  of 
positive  precursor  waves,  which  are  not  included  in  the  present  model.  Comparison  with  the  measurements  is,  thus, 
complicated  by  these  considerations  and  requires  further  study. 
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Figure  10.  Results  of  MSS  Shot  4286 
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Figure  10.  Results  of  MSS  Shot  4286  (Continued) 
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Figure  10.  Results  for  MSS  Shot  4286  (Continued) 


In  Figures  lOg  and  lOh,  pressure  gages  are  located  6  feet  above  the  bottom  interface.  In  these  cases  the  shapes  of 
the  wave  forms  differ  markedly  from  those  associated  with  the  2-inch  gage  positions.  The  figures  show  that  the  model 
is  in  qualitative  agreement  with  the  measurements,  but  the  pressure  wave  forms  and  impulses  in  these  cases  do  not  match 
the  measurements  as  well  as  when  the  gages  are  near  the  bottom.  In  Figure  1  Oh  an  unphysical  O  value  more  negative  than 
- 1  would  have  been  required  to  match  the  measured  impulse-time  curve.  Nevertheless,  the  results  demonstrate  that  the 
basic  mechanism  of  pulse  formation  is  probably  being  captured  by  the  model.  It  is  believed  that  the  agreement  could  be 
improved  with  model  refinement.  Because  oin  primary  interest  in  this  work  has  been  in  attacking  mines  at  locations  near 
the  bottom,  such  improvements  have  bear  deferred  imtil  a  later  time. 

Figure  1 1  shows  comparisons  betweai  calculations  and  measurements  for  the  discrete  Sabre  (D-Sabre)  line  charge 
system  in  1 0  feet  of  water  and  with  die  gage  positions  2  inches  above  the  bottom  interface.  Again,  calculations  are  shown 
in  red  and  measurements  in  blue.  It  should  be  noted  that  the  records  in  Figure  1 1  extend  over  longer  time  intervals  than 
those  shown  in  Figures  10a  through  lOf.  The  problems  with  late  time  impulse  match-up  evident  in  Figure  1 1  are  also 
characteristic  of  the  10-foot  M58  calculations,  but  these  times  are  not  displayed  in  Figures  10a  through  lOf. 
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Figure  11.  Results  for  D-Sabre  Shot  4287 
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Figure  11.  Results  of  D^Sabre  Shot  4287  (Continued) 
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The  D-Sabre  system  has  larger  charges  and  larger  spacings  between  charges  than  the  M58.  Apparently  because  of 
these  design  differences,  the  D-Sabre  pressure-time  history  differs  in  several  respects  from  the  M58  signature  (and  also 
from  that  of  continuous  Sabre,  which  is  closer  to  the  M58,  as  will  be  shown  below).  D-Sabre  has  fewer  bottom-reflected 
arrivals  between  the  leading  edge  of  the  pulse  and  surface  cutoff  time.  In  addition,  the  measured  pressures  show  broadened 
pressure  peaks,  distinctly  different  from  the  calculated  pressure  spikes,  that  are  possibly  a  result  of  finite  amplitude  effects. 
Because  fliere  fewer  arrivals  within  a  unit  interval  of  time,  differences  between  calculated  and  measured  arrival  times  are 
more  evident  for  tire  D-Saber  cases.  The  arrival  times  are  quite  sensitive  to  errors  in  gage  positions,  which  may  account 
for  some  of  the  observed  differences.  Anodier  source  or  significant  error  is  probably  the  assumption  of  sonic  wave  speeds. 

The  figures  verify  that  it  is  possible  to  select  constant  values  so  that  the  calculated  impulses  are  in  good 
agreement  with  the  measured  impulse  time  curves  up  to  the  time  of  surface  cutoff,  but  thereafter  the  calculations  fall  below 
the  measured  levels  which  rise  at  a  faster  rate.  The  time  of  surface  cutoff  is  not  nearly  as  distinct  in  the  D-Sabre 
measurements  as  it  appears  in  the  M5  8  records.  It  is  felt  that  this  may  be  associated  with  the  peak  broadening  mechanism 
and,  as  pointed  out  by  Britt,  with  the  arrival  of  doubly  reflected  (source-to-bottom-to-surface-to-gage  or  source-to- 
surface-to-bottom-to-gage)  signals  after  the  time  of  surface  cutoff  that  are  not  included  in  the  present  model.  The  absence 
of  multiply  reflected  arrivals  may  also  explain  the  poorer  late  time  agreement. 

Figure  12,  pertaining  to  the  continuous  Sabre  system,  completes  the  comparison  of  the  constant  O  model 
calculations  with  the  10-foot  water  depth  measurements.  Because  of  the  similarities  of  design  between  the  C-Sabre  and 
M58  S5^tems,  the  pressure  signatures  shown  in  Figures  12  resemble  those  of  Figures  10a  through  1  Of  The  problems  with 
late-time  match-up  are  again  apparent. 

One  concludes  from  all  10-foot  water  depth  comparisons  (i.e.,  from  Figures  10, 1 1,  and  12)  that  the  constant  ® 
model  provides  quite  reasonable  agreement  with  measured  inqjulse-time  curves  for  a  period  lasting  about  2  ps  after  arrival 
of  the  pulse  train.  It  must  be  remembered,  of  course,  that  O  values  have  been  selected  to  promote  this  agreement. 
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Figure  12.  Results  for  C-Sabre  Shot  4288 
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Figure  12.  Results  for  C-Sabre  Shot  4288  (Continued) 
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Figure  12.  Results  for  C-Sabre  Shot  4288  (Continued) 


Comparison  of  Constant  (D  Model  Calculations  with  Measurements  in  Water  3  Feet  Deep: 

Figures  i3,  14,  and  15  show  comparisons  between  calculated  and  measured  pressures  and  impulses  for  the  M58 
(shot  4283),  D-Sabre  (shot  4284),  and  C-Sabre  (shot  4285)  line  charge  systems,  respectively,  in  3  feet  of  water.  In  these 
shallow  water  shots  the  gages  were  positioned  closer  to  the  line  charge  than  in  the  10-feet-deep  cases.  Gage  A  was  at  a 
horizontal  standoff  of  6  feet  from  the  line  charge  axis,  gage  B  was  at  8  feet,  gages  C  and  D  were  at  10  feet,  gage  E  was 
at  13  feet,  gages  F  and  G  were  at  18  feet,  gages  H  and  I  were  at  23  feet,  and  gage  J  was  at  28  feet.  With  the  following 
exceptions  the  gages  were  located  2  inches  from  the  sand  bottom  in  the  plane  perpendicular  to  and  bisecting  the  line  charge 
a.xis.  The  D,  G,  and  1  gages  were  located  in  this  plane  at  the  indicated  standoffs,  but  at  a  height  of  2  feet  from  the  sand 
bottom  (one  foot  from  the  free  surface).  It  is  seen  that  the  F  gage  in  the  3-foot  water  depth  shots  was  at  the  same  position 
as  the  C  gage  in  the  10-foot  water  depth  shots,  but  that  all  other  gage  positions  were  different. 
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We  note  that  the  only  systematic  difference  between  M58  shots  4283  and  4286,  between  D-Sabre  shots  4284 
and  4287,  and  between  C-Sabre  shots  4285  and  4288  is  the  water  depth  (see,  e.g..  Table  I).  This  means  that  up  until  the 
arrivals  of  the  first  fi-ee  surface  reflections,  these  paired  shots  can  be  regarded  as  identical  except  for  random  differences. 
The  pre-surface-cutoff  portions  of  the  histories  shown  for  the  3-foot  depth  cases  can  be  thought  of  as  replicate 
observations  (in  most  cases  fi-om  different  gage  positions)  of  the  pressure  fields  generated  in  the  10-foot  water  depth 
shots.  We  will  use  this  fact  in  the  next  section  to  make  a  comment  on  the  reproducibility  of  the  Eglin  tests.  It  is  also  useful 
to  keep  this  fact  in  mind  while  examining  Figure  1 3,  where  the  early  time  M5  8  pressure  histories  are  shown  on  expanded 
scales  and  the  structures  of  the  leading  spikes  are  more  apparent  than  before. 

In  Figure  13  (the  M58  linear  charge  array)  surface  cutoff  occurs  approximately  0.25  to  0.45  ps  after  the  first 
arrival  of  the  pulse  train  depending  upon  gage  position.  This  is  contrasted  with  approximately  1.2  to  2.2  ps  before  surface 
cutoff  for  the  10-foot  water  depth  cases.  Aside  fi'om  the  initial  spike,  the  buildup  of  the  pressure  ramp  appears  to  be 
reasonably  well  approximated  by  the  constant  O  model  calculations.  The  calculated  impulse-time  curve  is  also  in  apparent 
agreement  with  the  measurements  over  the  time  span  shown  in  the  figures.  In  fact,  however,  the  impulse  falls  below  the 
measurements  at  later  times  for  the  M58  cases  and  also,  as  will  be  seen,  for  the  shallow  continuous  Sabre  cases. 
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b.  Gage  B,  8-Foot  Standoff,  2-Inch  Height,  <t>  =  -0.982 


Figure  13.  Results  for  MSS  Shot  4283 
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c.  Gage  C,  10-Foot  Standoff,  2-Inch  Height,  0  =  -0.985 
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d.  Gage  E,  13-Foot  Standoff,  2-Inch  Height,  0  =  -0.99 


Figure  13.  Results  for  MSS  Shot  4283  (Continued) 


Figures  14  (the  D-Sabre  linear  array)  and  15  (the  C-Sabre  linear  array)  demonstrate  the  behavior  of  the  model 
out  to  later  times.  At  present  this  behavior  is  not  fully  imderstood.  The  D-Sabre  impulse  calculations  are  in  excellent 
agreement  with  the  experimental  measurements,  although  the  agreements  between  the  pressure-time  curves,  from  which 
the  impulse-time  curves  were  integrated,  are  less  convincing.  Conversely,  the  C-Sabre  (and,  although  not  shown,  the  M58) 
calculated  impulses  lie  considerably  below  the  measurements  at  later  times,  while  the  pressiue-time  curves  at  early  times 
are  in  quite  reasonable  agreement.  To  obtain  these  results,  calculations  were  stopped  after  the  arrival  of  the  20th  free- 
surface-reflected  wave  in  Figure  14  (and  also  Figure  13)  and  after  the  arrival  of  the  120th  free-siuface-reflected  wave  in 
Figure  15.  Six  times  the  number  of  free-surface-reflected  arrivals  were  necessary  in  the  (Figure  15)  C-Sabre  calculations 
to  create  record  lengths  comparable  to  the  D-Sabre  calculations  because  the  C-Sabre  array  contains  six  times  the  number 
of  charges  as  contained  in  the  D-Sabre  array.  The  agreement  between  the  C-Sabre  calculated  and  measured  impulses 
seems  to  improve  as  the  standoff  from  the  line  charge  axis  is  increased.  While  the  underprediction  of  late  time  impulse 
by  the  M58  and  C-Sabre  models  might  be  explained  by  the  absence  of  multiple  reflection  arrivals  as  suggested  by  Britt, 
it  is  not  understood  why  the  D-Sabre  model  performs  so  well  without  the  multiple  reflections  included. 
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f.  Gage  H,  23-Foot  Standoff,  2-Inch  Height,  0  =  -0.985 


Figure  14.  Results  of  D-Sabre  Shot  4284  (Continued) 


24 


PRESSURE  (PSi) 


IHTR2058 


a.  Gage  A,  6-Foot  Standoff,  2-Inch  Height,  d)  =  -0.96 
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c.  Gage  C,  10-Foot  Standoff,  2-Inch  Height,  0  =  -0.985 


Figure  15.  Results  for  C-Sabre  Shot  4285 
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d.  Gage  E,  13-Foot  Standoff,  2-Inch  Height,  (J> «  -0.99 
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f.  Gage  J,  28-Foot  Standoff,  2-Inch  Height,  O  =  -0.999 

Figure  15.  Results  for  C-Sabre  Shot  4285  (Continued) 
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MODEL  CONSISTENCY  AND  MEASUREMENT  REPRODUCIBILITY 


It  was  noted  earlier  that  pressure  time  curves  for  3-  and  10-foot  water  depth  cases  of  the  same  linear  charge  array 
system  should  be  identical  for  gages  at  the  same  positions  before  the  arrival  of  surface-reflected  waves.  This  fact  provides 
a  means  of  judging  the  reproducibility  and  consistency  of  the  measurements  made  in  the  1995  Eglin  test  series  and, 
unexpectedly,  also  of  the  modeling  analysis.  Figure  16  shows  the  records  and  calculations  associated  with  the  F  gage  in 
D-Sabre  shot  4284  and  the  C  gage  in  D-Sabre  shot  4287.  Shot  4284  was  at  a  water  depth  of  3  feet  and  shot  4287  was 
at  10  feet.  Both  gages  were  at  a  standoff  of  18  feet  and  were  2  inches  above  the  sand.  Unfortunately,  because  of  the 
presence  of  hi^  frequency  electrical  interference  in  the  F  gage  records  of  the  other  shallow  shots,  these  are  the  only  two 
records  that  can  be  so  compared.  No  other  gages  were  placed  at  identical  positions. 


0.006  0.0062  0.0064 
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Figure  16.  Comparison  of  identicai  Gage  Position  Resuits 
for  D-Sabre  Shots  4284  (Left)  and  4287  (Right) 


Figure  16  reveals  quite  close  agreemoat  between  the  two  experimentaUy  measured  pressure  histories  at  early  times 
(shown  in  blue).  The  second  peaks  are  particularly  comparable.  The  disagreement  between  leading  peak  pressures  (and 
also  of  the  second  peak  pressures)  may  be  a  result  of  the  1-ps  gage  sampling  rate  used,  reflected  in  the  visible  stepwise 
fine  structures,  which  appears  to  be  marginally  adequate  to  effect  resolution  of  the  peaks.  The  absence  of  the  third  peak 
in  the  left  figure  can  be  attributed  to  surface  cutoff,  which  always  occurred  in  advance  of  the  calculated  cutoff  time  due 
to  the  supersonic  speeds  of  the  actual  waves.  Based  on  the  agreement  of  the  pressures  of  the  second  peaks,  one  might 
tentatively  conclude  that  pressures  from  linear  charge  arrays  at  Eglin  would  be  reproducible  from  shot  to  shot.  The  two 
D-Sabre  shots  were  conducted  one  week  apart  and  in  separately  prepared  sites  of  the  Eglin  test  area. 

The  consistency  of  the  modeling  analyses  apphed  to  both  cases  is  seen  in  Figures  11c  and  14e,  where  it  is  found 
that  the  reflection  factors  selected  by  matching  impulse  time  curves  were  both  $  =  -0.978,  i.e.,  identical  to  the  third 
decimal  place.  These  assignments  were  earned  out  independent^  and  before  it  was  realized  that  the  comparisons  in  Figure 
16  could  be  made.  The  fact  that  the  impulse-time  curve  in  Figure  1  Ic  is  quite  dissimilar  in  shape  from  that  in  Figure  14e 
due  to  the  difference  in  surface  cutoff  times  lends  confidence  to  the  modeling  approach. 
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Because  the  same  $  values  were  used,  the  calculated  pressures  in  both  shots  were  identical  at  early  times. 
Following  surface  cutoff,  the  pressure  ramp  in  the  left  graph  of  Figure  16  is  eliminated  by  the  arrival  of  surface- 
reflected  tension  waves.  The  pressure  spikes  associated  with  the  direct  and  bottom-reflected  arrivals,  however,  remain. 
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DEVELOPMENT  OF  THE  VARIABLE  0  MODEL 


The  strategy  followed  for  development  of  a  variable  $  model  was  to  associate  each  constant  model  reflection 
coefficient  value  appearing  in  Figures  10  through  15  with  an  appropriate  reflection  angle  y  and  determine  functional  forms 
consistent  with  all  such  fitted  points  that  could  be  used  in  file  variable  $  model  code.  In  the  code  different  reflection 
coefficient  values  would  then  be  assigned  to  the  bottom-reflected  ray  of  each  charge  of  the  system  by  inserting  the 
reflection  angle  of  the  into  the  fimctional  form.  The  comphcation  in  this  approach  is  that  each  impulse  time  curve  and 
fitted  O  value  is,  in  fact,  associated  with  a  range  of  reflection  angles.  Nevertheless,  this  information  was  useful  in  choosing 
final  functional  forms  that  seemed  to  produce  suitable  overall  matches  with  the  data.  While  a  more  formal  statistical 
procedure  for  estimating  the  functional  forms  could  certainly  have  been  pursued,  the  approach  taken  seenied  more  direct 
and  less  demanding  fi-om  a  programming  standpoint  with  the  Mathcad  (Version  6+)  software  being  used. 

An  examination  of  the  reflection  coefficient  values  (<&)  and  reflection  angle  (y)  ranges  for  all  six  shots  suggested 
fliat  eadi  line  charge  system  required  a  separate  functional  relationship  between  O  and  y  to  produce  the  best  agreement 
with  the  data.  The  underlying  reason  for  this  is  not  completely  clear.  But  because  the  levels  of  sand  compaction  caused 
by  the  explosions  (which  could  affect  the  values  of  O)  would  be  different  for  each  charge  system  due  to  differences  in 
explosive  type,  charge  weight,  and  system  configuration,  the  use  of  separate  fimctional  forms  appears  reasonable. 

The  fitted  leflecfion  coefficiaits  and  reflection  angle  ranges  associated  with  the  M58,  D-Sabre,  and  C-Sabre  systems 
are  shown  in  Figures  17, 18,  and  19  respectively.  In  each  of  these  the  y  ranges  appear  as  horizontal  lines  bound^  by  o 
symbols  if  the  measurement  was  associated  with  a  10-foot  water  depth  case  and  x  symbols  if  it  was  associated  with  a  3- 
foot  water  depth  case.  In  the  latter  the  horizontal  lines  are  also  heavier  in  weight  so  that  they  can  be  easily  distinguished 
\^hen  overlapping  with  the  10-foot  case  lines  occurs.  The  end  points  of  each  range  represent  the  maximum  and  minimum 
values  of  y  for  bottom-reflected  waves  that  arrive  between  the  leading  edge  of  the  pulse  train  and  the  time  of  surface 
cutoff.  This  measure  of  range  was  chosen  because  it  reflected  the  portion  of  the  pressure  history  best  fitted  by  the  model. 

The  blade  curves  appearing  in  Figures  17, 18,  and  19  are  quadratic  forms,  suggested  by  the  data,  that  were  observed 
to  provide  the  best  overall  fits  to  the  impulse  time  data  using  tihe  variable  $  model.  The  specific  equations  used  in  each 
case  are  indicated  below; 


^M58  =  14y2  +  Oy  -  1  (16) 

®D-Sabre  =  22y2  +  0.2y  -  1  (17) 

%sab«  =  15y2  +  O.ly  -  1  (18) 
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Appendix  A  contains  the  results  of  calculations  made  with  the  variable  ®  model  in  which  the  reflection  coefficients 
were  determined  for  each  individual  charge  and  gage  position  using  Equations  16, 17,and  18.  Comparison  of  Appendix 
A  with  Figures  10  through  15  shows  that  the  variable  $  model  sli^tly  underperforms  the  constant  model  in  terms  of 
overall  agreement  with  the  measured  impulses.  But  in  most  cases  the  performances  are  comparable. 

Figures  17, 18,  and  19  indicate  that  the  best  agreement  of  the  variable  ^  model  pressures  with  the  measured  data 
was  obtained  by  associating  the  fitted  $  values  (determined  with  the  constant  $  model)  with  the  highest  values  of 
reflection  angles.  That  is,  the  best  fit  functional  forms  are  located  along  the  right  edges  of  the  plotted  ranges  of  y  in  all 
three  figures.  The  initial  expectation  was  that  the  best  fit  curves  would  pass  through  mi^ange  values  that  were 
representative,  in  some  average  sense,  of  the  overall  fit.  It  is  suggested  that  the  largest  values  of  y  provide  the  best  fits 
b^use  they  are  the  values  associated  with  the  earliest  and  most  rapid  rise  of  the  pressure  ramp  and  to  which  the  overall 
shape  of  the  ramp  is  most  sensitive. 


Figure  17.  MSS  Reflection  Coefficients  Versus  Angie 
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Figure  18.  D-Sabre  Reflection  Coefficients  Versus  Angle 
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REFLECTION  ANGLE  (RADIANS) 

Figure  19.  C-Sabre  Reflection  Coefficients  Versus  Angle 
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CONCLUSIONS 


It  was  found  that  ray  tracing  modeling  techniques  could  be  effectively  used  to  represent  the  pressure  histories 
generated  by  linear  arrays  of  explosive  charges.  It  is  concluded  that — 

•  Ray  tracing  techniques  can  be  used  to  provide  rapid  computations  of  explosive  array  pressure  histories 
wiA  accuracies  that  are  sufficient  for  the  analysis  of  mine  structural  responses. 

•  It  is  expected  that  the  method  would  also  be  effective  for  representing  the  pressures  associated  with  odier 
charge  configurations  and  geometries  such  as  DET  arrays,  nonlinear  arrays,  and  arrays  on  sloped 
bottoms. 

•  The  theory  of  tiiis  modeling  approach  should  be  pursued  to  permit  the  reflection  coefficient  to  be 
expressed  as  a  function  of  water  and  sand  mechanical  properties. 
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Appendix  A 

COMPARISON  OF  VARIABLE  <t>  MODEL  CALCULATIONS  WITH  MEASUREMENTS 
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COMPARISON  OF  VARIABLE  O  MODEL  CALCULATIONS  WITH  MEASUREMENTS 


Figures  A-1  through  A-6  show  comparisons  of  calculations  made  using  the  variable  ®  model  with  measur^ents 
for  the  conditicms  of  the  1995  Egjin  mine  vulnerability  tests.  ^  The  numeral  of  each  figure  designation  pertains  to  the  shot 
order  followed  in  die  Eglin  tests.  As  an  exan^le.  Figures  A- la  through  A- If  show  measurements  recorded  during  the  first 
shot  of  the  test  series.  The  ending  letter  in  the  designation  denotes  the  cable  and  gage  as  explained  in  the  text.  Blue  curves 
are  the  measurements  and  red  traces  are  the  calculations.  The  figure  titles  indicate  the  line  charge  array  system,  the  shot 
number,  and  the  gage  position. 
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b.  Gage  B,  8-Foot  Standoff,  2-Inch  Height 

Figure  A-1.  Results  for  MSS  Shot  4283  with  Variable  (p 
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f.  Gage  J,  28-Foot  Standoff,  2-Inch  Height 

Figure  A-1.  Results  for  M58  Shot  4283  with  Variable  0  (Continued) 


a.  Gage  A,  6-Foot  Standoff,  2-Inch  Height 


b.  Gage  B,  8-Foot  Standoff,  2-Inch  Height 

Figure  A-2.  Resuits  for  DSabre  Shot  4284  with  Variabie  0 
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c.  Gage  C,  10-Foot  Standoff,  2-Inch  Height 


d.  Gage  E,  13-Foot  Standoff,  2-Inch  Height 


e.  Gage  H,  23-Foot  Standoff,  2-Inch  Height 


Figure  A-2.  Results  for  D-Sabre  Shot  4284  with  Variable  0  (Continued) 
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a.  Gage  A,  12-Foot  Standoff,  2-lnch  Height 


b.  Gage  B,  15-Foot  Standoff,  2>lnch  Height 
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c.  Gage  C,  18-Foot  Standoff,  2-Inch  Height 


Figure  A-4.  Results  for  MSS  Shot  4286  with  Variable  0 
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d.  Gage  E,  22-Foot  Standoff,  2-lnch  Height 
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e.  Gage  H,  30-Foot  Standoff,  2-Inch  Height 
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f.  Gage  J,  35-Foot  Standoff,  2-Inch  Height 


Figure  A-4.  Results  for  MSS  Shot  4286  with  Variable  0  (Continued) 
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b.  Gage  B,  15-Foot  Standoff,  2-Inch  Height 
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c.  Gage  C,  18-Foot  Standoff,  2-Inch  Height 


Figure  A-5.  Results  for  D-Sabre  Shot  4287  with  Variable  0 
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d.  Gage  E,  22-Foot  Standoff,  24nch  Height 


e.  Gage  H,  30-Foot  Standoff,  2-)nch  Height 


f.  Gage  J,  35-Foot  Standoff,  2-inch  Height 

Figure  AS.  Results  for  D-Sabre  Shot  4287  with  Variable  0  (Continued) 
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a.  Gage  A,  12-Foot  Standoff,  2-Inch  Height 


b.  Gage  B,  15-Foot  Standoff,  2-lnch  Height 
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c.  Gage  C,  18-Foot  Standoff,  2-lnch  Height 


Figure  A-6,  Results  for  C-Sabre  Shot  4288  with  Variable  0 
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d.  Gage  E,  22-Foot  Standoff,  2-Inch  Height 
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e.  Gage  H,  30-Foot  Standoff,  2-Inch  Height 


f.  Gage  J,  35-Foot  Standoff,  2-Inch  Height 


Figure  A-6.  Results  for  C-Sabre  Shot  4288  with  Variable  0  (Continued) 
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